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Abstract 

The existence of an anomalous U(l) symmetry is shown to play a crucial role in the 
supersymmetric radiative seesaw model for neutrino masses. It explains the smallness of 
some couplings related to neutrino mass generation in a favorable way in addition to cause 
the hierarchical structure of Yukawa couplings of quarks and leptons. If it is spontaneously 
broken to a Z2 subgroup, this Z2 symmetry can make a lifetime of the lightest field with 
its odd parity extremely long. Thus, the model has an additional dark matter candidate 
other than the lightest neutralino, which appears in the i?-parity conserved MSSM. We 
discuss the nature of dark matter by taking account of its relation to the neutrino mass 
generation and the lepton flavor violation. 
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1 Introduction 



The recent astrophysical observations pQ and neutrino oscillation experiments [2] require 
that the standard model (SM) should be extended so as to include dark matter (DM) 
[3] and small neutrino masses. A radiative seesaw model proposed by Ma [1] is a simple 
and interesting possibility among such extensionsjf] In this model the existence of DM is 
intimately related to the neutrino mass generation. The combined study of the neutrino 
oscillation data, the lepton flavor violating processes such as \i — > and the DM relic 
abundance can give strong constraints on the model [5J E|. On the other hand, the recently 
reported cosmic ray anomaly [12j IT3] attracts a lot of attention as a target of the DM 
physics. Various work relevant to this has been done on the basis of the annihilation of 
DM [H [13 d3 and also the decay of DM US SB 123 EH |22J. If we impose that 
this neutrino mass model should explain this anomaly, it could bring the model additional 
valuable information not only on the mass and the interaction of DM but also on the lepton 
flavor structure [17J. Although the model shows interesting and promising features on 
phenomenology, the model cannot give any answer to the hierarchy problem unfortunately. 
Previous models are usually constructed in the nonsupersymmetric framework. Thus, the 
detailed study on the supersymmetrization of the model is one of the remaining subjects 
relevant to the model. 

The supersymmetric extension of the model has been considered in Two 
DM candidates appear in these extensions as long as i?-parity conservation is assumed. 
The existence of two kinds of Z 2 symmetry, that is, i?-parity and Z 2 which forbids tree- 
level neutrino masses, guarantees the stability of these two DM candidates. In that 
case both of them contribute to the expected DM relic abundance^ However, if one of 
them is unstable but has an longer life time than the age of universe, this decaying DM 
may cause the additional contribution cosmic ray, which brings anomalous excess of the 
charged particle flux over the expected background. Such a possibility has been studied 
in a supersymmetrized model in [21], where Z 2 is assumed to be weakly broken by the 
anomaly effect. Since the anomaly induced interaction is strongly suppressed, it causes 

3 There are a lot of study for the DM nature in radiative neutrino mass generation models and other 

phenomenological features of such models [51 [71 [HI HI EH HI] • 

4 A relevant supersymmetric model is considered in a different context in [25) . 
5 Multi-component DM models have been considered in several contexts [2"B]. 
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the DM decay but its life time can be long enough. 

In this paper, we consider a modification of the original one given in [21] by intro- 
ducing the an anomalous U(l) gauge symmetry. Spontaneous breaking of this symmetry 
causes the weakly broken Z 2 symmetry at low energy regions naturally in addition to the 
conserved .R-parity. It makes us possible to understand the smallness of some couplings 
through the Froggatt-Nielsen mechanism [271 128] , which is required to generate the small 
neutrino massesjj Simultaneously, hierarchical masses of quarks and charged leptons are 
also explained by the same origin. The model also shows another interesting features 
relevant to the DM phenomenology. We give its analysis which includes the cosmic ray 
anomaly expected in this model, the characteristic gamma predicted in the DM radiative 
decay and the direct search of the stable DM through the elastic scattering with nuclei. 

The paper is organized as follows. In section 2 we give the low energy effective model 
with an anomalous U(l) gauge symmetry. Spontaneous breaking of this symmetry gives 
the Z 2 symmetry which is relevant to the neutrino mass generation. The neutrino mass 
and mixing, the lepton flavor violating processes and the DM relic abundance are studied 
in this framework. In section 3 the nature of DM is discussed from a view point of both 
direct and indirect searches. Section 4 is devoted to the summary. 



2 A supersymmetric extension 

The radiative seesaw model proposed in [1] is an extension of the SM with three right- 
handed neutrinos and an inert doublet scalarJl) It has no coupling with quarks and no 
vacuum expectation value (VEV). Although the model has interesting phenomenological 
features as discussed in El [8], it can not give an answer to the hierarchy problem. 
Thus, it is worthy to consider the supersymmetric extension of the model and to study 
phenomenology associated to it. The supersymmetrization of the model requires to in- 
troduce inert doublet chiral superfields rj u and r\ d and also a singlet chiral superfield to 
the minimal supersymmetric SM (MSSM) [21]. They play the similar role as the inert 
doublet scalar in the original nonsupersymmetric model and bring required terms in the 

6 An anomalous U(l) symmetry has been considered in various phenomenological contexts |29j . In a 
radiative seesaw mechanism, it is also considered in [30j . However, the present model is different from it 

in the nature of dark matter and also the i?-parity conservation. 

7 The origin of the second doublet is discussed in [31] , 
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Table 1: Matter contents and their quantum number. X represents the charge of the anomalous X5(l)x 
and each n<sr is an integer. Z 2 is a remnant symmetry of U(l)x caused by the symmetry breaking due to 
<£±}^0. 



Lagrangian for the neutrino mass generation. As the origin of the Z 2 symmetry which 
forbids tree level neutrino masses, we suppose the existence of an anomalous U(l)x gauge 
symmetry. We assume that the .R-parity is conserved. Matter contents of the model and 
their relevant quantum numbers are summarized in Table 1. 

2.1 A low energy effective model 

First, we discuss spontaneous breaking of the anomalous U(l)x symmetry at a high energy 
scale, which brings the low energy effective theory. Hierarchical couplings and masses are 
found to be generated in the low energy theory through this breaking. The vacuum is 
expected to be determined as the flat direction of the .D-term of U(l)x- The .D-term for 
relevant hidden fields is given by 
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£x is the U(l)x Fayet-Iliopoulos .D-term. It is expressed in string models as 

TV Y" 

where M p \ is the reduced Planck mass and X stands for the anomalous U(l)x charge 
of the fields. On the other hand, the lowest order superpotential for these fields can be 
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written as 

W h = -^-E + E n+ + -^0 2 E 2n ^ +1 , (3) 

Mir + - m p 2 ^ " 

where c + and are considered as real constants. If we suppose supersymmetry breaking 
in the hidden sector, supersymmetry breaking terms appear in the scalar potential of the 
hidden sector. The scalar potential derived from Wh may be written as 

C|(2n, + 1) 2 2w / ^ +C »(2^+1) ! 2 2n 

M 4„, ~ ^ M «++2n,-2 v + - ; r - 

-FxS 2 ! - FsS 2 , , (4) 

where F 1 and F 2 represent the VEVs of F-components of some chiral superfields in the 
hidden sector. They bring the supersymmetry breaking of O(10 n ) GeV in the hidden 
sector and also induce the soft terms of 0(1) TeV through the gravity mediation in the 
observable sector. 

Here we assume that the VEVs of E± and are real, for simplicity. Then, the min- 
imization of the potential V — V D + V F along the D-flat direction suggests that there 
exists a vacuum defined by 

(E_)^^>><E + >, (0) = O, (5) 

as long as c + is sufficiently suppressed. For example, if TrX ~ 150/ g\ is satisfied, 
(E_) ~ 0.2M p i is expected. Moreover, (E + ) ~ 10" 4 M p i is also expected for sufficiently 
suppressed values of c + such as O(10~ 7 ). Although this kind of vacuum can be realized 
only for the finely tuned parameters, we assume it in the following study. 

If we adopt this vacuum to fix the low energy effective theory, the superpotential 
invariant under the imposed symmetry is obtained with the effectively induced parameters 
such as 

W = hYjQiUjH^ + h^QiDjHd + hfLiEiHd + fiHH^, 

+ hfj LiNji] u + \ u r] u H d (j) + \ d r] d H u (j) + /inVuVd + ^M; AyVj + -^(p 2 . (6) 

The invariance of each term under _RxU(l)x is guaranteed since the effective couplings 
and masses are generated through the VEVs (E±) as shown below (see also Appendix 



A). The MSSM superpotential is contained in the first line, while the second line includes 
additional terms to the MSSM. Yukawa couplings for the charged leptons and the mass 
matrix for the right-handed neutrinos are supposed to be flavor diagonal jf| 

The effective parameters in this superpotential are induced from higher order invariant 
interaction terms which contain a suitable number of E±. The dominant contribution are 
determined by the lowest order term of the following form: 



The original coupling constants and in the nonrenormalizable interaction terms 
are considered to be values of 0(1). If the singlet scalars E± obtain the VEVs as discussed 
above, these VEVs cause hierarchical structure in the Yukawa couplings of quarks and 
leptons, and also suppress several parameters in the superpotential as found from eqs. (JSJ) 
and (j7J). In fact, as long as suitable U(l)x charges are assigned to quarks and leptons, we 
find that the hierarchical mass eigenvalues and mixing are generated via Froggatt-Nielsen 
mechanism [27]. In the similar way, the Yukawa couplings are largely suppressed, 
and also //,, and Mj can take values of 0(1) TeV. These parameters are relevant to the 
neutrino mass generation and their values can be favorable for it as seen in the next 
part. We give more detailed discussion by giving such concrete examples for the charge 
assignment in the Appendix A. We also find that (E±) breaks the anomalous U(l)x 
symmetry spontaneously to its Z2 subgroup. This is obvious from the charge assignment 
shown in Table 1. This Z2 symmetry plays a crucial role in the DM phenomenology. 

Here we also note an interesting point related to the anomaly induced interaction. 
We have as an only renormalizable operator which breaks the XJ(l)x invar iance but 
is invariant under the SM gauge symmetry and the i?-parity. It is not included in the 
low energy superpotential ([6]) since we can not make it XJ(l)x invariant by multiplying 
any number of E±. However, it can be U(l)x invariant if it is accompanied by a dilaton 
chiral superfield S. As is well known, if the anomalous U(l)x gauge transformation 
Vx — » Vx + i(A — At)/2 is associated with the shift of the dilaton field such as 



This can be justified as long as the relevant original couplings and j/y in eq. ([7]) are flavor 




(7) 



S ^ S + id GS A 
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diagonal. 
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the anomaly cancellation for XJ(l)x is completed [33| IM] . Noting this fact, we observe 
that the following non-perturbative superpotential is also invariant under the imposed 
symmetry [31] 1^1 

W np = CiM pl e- bl Lir] u , 

(2nf + 2n Vu + l) G 1927r 2 (2nf + 2ry + 1) 

6i = 5 (9) 

where we use (S) and q are constants of 0(1). The expression for hi shows that it is 
determined by the U(l)x charge of the field contents of the model including the hidden 
sector. We also note that this term breaks the remnant Z 2 symmetry very weakly. The 
scale of its violation is determined by both the value of TrX and the charges of Lj and 

Vu- 

Soft supersymmetry breaking terms associated with the superpotential W and W np 
are introduced as follows, 

+A(hf j L i N j f) u + \ u fj u H d (f> + \ d rj d H u 4> + h.c.) 

-B (^ H H u H d + f^fjjja + + ^MiN? + aM&e-^Lifjv + h.c. J , (10) 

where the additional part to the MSSM is listed alone. The scalar component is repre- 
sented by putting a tilde on the character of the corresponding chiral superfield except for 
the ordinary Higgs chiral superfields H u and H d . Universality of the soft supersymmetry 
breaking A- and 5-parameters is assumed to be satisfied here. Moreover, we confine our 
consideration to the case in which the soft scalar masses for all scalar partners of quarks 
and leptons are flavor diagonal and universal. They are denoted by mo in the following 
analysis and we assume A — B — mo, for simplicity. 

We should remind the reader that the soft scalar masses have the contributions from 
the anomalous U(l)x D-term such as 

Am 2 = &X t [C h + 2n Hu (H° u ) 2 + 2n Hd (H° d ) 2 ) , (11) 

where Ch stands for the remnant contribution from the hidden sector. Although these are 
flavor dependent, we confine our study to the case m^ S> Am- , which may be realized on 



'The effect of the appearance of this term has been studied in the case of discrete symmetry in [24] . 
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the vacuum with a finely tuned Ch- In such a case there is no new dangerous origin for the 
flavor changing neutral current (FCNC) caused by the soft supersymmetry breaking terms. 
In the following study the effective parameters in eq. ([6]) and also the soft supersymmetry 
breaking parameters in eq. (fTUl) are treated as real except for hfj, for simplicity. 

2.2 Neutrino mass and mixing 

Neutrino masses are generated through the one-loop diagram shown in Fig. 1 when the 
Higgs doublet scalars H® and H% obtain the VEVs. As found from eq. (j4"0l) . since X u and 
Xd are very small, the mixing between (fju*,Vd) an d can be treated as an insertion 

in the calculation of these diagrams with good accuracy. The mass matrices for (??°*,?7°) 
and (<t>*,(j)) are written as 

Al„„ = I < M , M. = U M * ^ I (12) 



B N ml ) 2 1 m 2 



where m 2 ^ d ~ m\ + /i 2 + \\ A v 2 d u and m| ~ m\ + /i| + \\v \ + X 2 d v 2 . If we define the mass 
eigenstates of these mass matrices by 



cos 9 V shifty \ I Vu* \ ( < f ) + i / cos Off, sin^ 






- sin 9 n cos 9 V I \ ri i I \ 4>- I \ — sin 9^ cos 6*^ 
the mass eigenvalues and the mixing angles can be written as 

<± = \ (< + < ± v « - O 2 + AB2 ^ ) ' tan2 ^ = ^r^' 

7T 

™5± = ™\ ± ^ = -. (14) 

We find that the scalar superpartners (N*, Ni) of Ni satisfy the same relations as the ones 
of (0*, <fi) as shown above. Their mass eigenvalues Mf ± can be read off from the expression 
for m 2 ± by replacing ^ and m| with Mj and j«q + M 2 , respectively. 

We can calculate the neutrino masses generated through the one-loop diagrams by 
using these. Since is expected to be larger than /i 2 and m§, we find that the dominant 
contribution is caused by the diagram (a) in Fig. 1 and it is estimated as 



(M v ) a p = ^ u ^ dV ^^ m 2 ^ haihpiMi (g{M h m v+ ) cos 2 9 V - g{M h m v _) sin 2 9 V 

i=i 

-f(M i ,m v+ ,m v -)cos29 11 ), (15) 



(a) 



(b) 



(c) 



+ 




Fig. 1 One-loop diagrams contributing to the neutrino masses. Two diagrams (b) and (c) have an 
internal line of the scalar 4> or its fermionic partner, respectively. 

where / and g are defined as 

m 2 a m 2 ln(m^/m„) + m 2 m 2 c ln{ml/m 2 ) + m 2 c m\ In {m 2 a /m 2 c ) 



f(m a} m b} m, 
g{m a ,m b ) = 



{ml - m 2 a ) {ml - m 2 ){m 2 - m 2 c ) 
m 2 — m 2 . + m\ \w{m 2 a j m\) 



m 



2\2 



(16) 



Other two diagrams (b) and (c) with the component of 4> as an internal line are expected 
to bring the subdominant contributions. We give the explicit expressions of their contribu- 
tions in the Appendix B. The universal soft supersymmetry breakings make the situation 
simple since 8 V = tt/A is satisfied. If Mi and m Vu d have the values of 0(1) TeV, the mass 
eigenvalues of neutrinos must be controlled by a small parameter X u Xd of 0(1CT 8 ). As 
long as the anomalous U(l)x charge is assigned suitably, these can be naturally realized 
as found from the examples given in the Appendix A. 

Here we focus our attention to the flavor structure of the neutrino Yukawa couplings^ 

Em 



h 



N 



0, h% = h 



,N 



.,„ Tl - hi a = i, 2), h» = h% = -h& = h. (i7) 

This flavor structure is consistent with the U(l)x invariance as long as all of take 
the same values except for y^, which should be assumed to be y^ — for i — 1, 2 (see 
Appendix A). An interesting point of this flavor structure is that the neutrino mass matrix 
in eq. (H5p takes the following simple form: 

/ \ / 1 



M, 



1 1 
\0 1 1/ 



(^Ai + h 2 2 A 2 ) + 



-l\ 

-1 

1 / 



h 2 3 A 3 , 
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This structure is chosen adhoc here since it has several interesting features as shown below. 
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where A,- fixes the mass scale for the neutrino masses as follows 11 



A,- 



Av 2 M 



g{Mi, m v+ ) - g(Mi, m n _) 



A 



A^A^tan/3 



(19) 



327r 2 V ^~n-^j ,-, r '-"i + tan 2 /? 

where sin (3 = v u /v and cos (3 = v^/v. This mass matrix automatically derives the tri- 
bimaximal mixing, which is favored by the neutrino oscillation data. In fact, it is easily 
checked that the MNS matrix for this neutrino mass model is given by 

\ / 1 \ 

u MNS = =h -4 -4 



\ 4= — — / 

\ V6 V3 s/2 / 

where Majorana phases a ii2 are expressed as 

1 
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(20) 
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2 A X cos 2ipi + 
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2 A 2 cos 2ip 2 



(21) 



by using ifi = arg(/ij). Here it should be reminded that Majorana phases a\ t2 do not 
affect the neutrino oscillations. 

Taking account of the nature discussed above and also the fact that one of the eigenval- 
ues of the mass matrix (fT8l) is zero, the remaining mass eigenvalues are found to be equal 



to a/ Am 2 tm and \/Am s 2 ol . This is required for the explanation of neutrino oscillation 
data. Thus, we find that the model parameters should satisfy the relations 
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hlk 



2 iy 2 \ 
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3^3| 



a/ Am 
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(22) 



The neutrino Yukawa couplings hi and the right-handed neutrino masses M, should cor- 
relate each other so as to satisfy these relations. Phenomenological study of the model 
should be proceeded under these constraints. In the following discussion, we restrict our 
study to the case with Mi ^ M 2 < M3, which allows us to take Ai ~ A 2 o Thus, free 
parameters relevant to the analysis of DM phenomenology are summarized as 



Mi, M 3 , A, fiy, m . 



(23) 



11 If the mass matrix is modified from eq. (|15j) due to other contributions neglected as the sub-dominant 
ones here, the ambiguity is confined into the Ai. The MNS matrix is not affected as long as the condition 
(fTT| is satisfied. 

12 It may be useful to note that this case has some advantages as discussed in [3 [17], but not only for 
simplicity. 
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Fig. 2 One-loop diagrams causing the lepton flavor violating processes £ a — > £^7. 



If we suppose fixed the values for Mi, fi v and mo, the conditions fl22|) determine a value 
of A a/ I h\ + ftjj and A | ^2.3 1 . A larger A gives a smaller value of neutrino Yukawa couplings 
\\J\h\ + h\\ and A | 1 . This feature becomes crucial when the constraints from lepton 
flavor violating processes and DM relic abundance are taken into account. 

2.3 Constraints from lepton flavor violating processes 

Since we suppose the slepton mass matrix is flavor diagonal and universal, there are no 
new FCNC source in the slepton sector. The FCNC is induced only through the MNS 
matrix elements which appear in the Higgsino exchange diagrams. Thus, if Higgsinos 
are much heavier than gauginos, the ordinary contributions to the FCNC caused by the 
supersymmetric partners are sufficiently suppressed in this model. On the other hand, 
the extension for the neutrino mass generation induces the one-loop contribution to the 
lepton flavor violating processes such as fi — > ey as in the nonsupersymmetric case. Their 
diagrams are shown in Fig. 2. Since these diagrams do not need the mixing between r\ u 
and r)d unlike the case of neutrino masses (see Fig. 1), it causes the large contributions 
to these processes. In that case the dominant contributions to the lepton flavor violating 
processes are given by these diagrams. They give the constraints on the model, which is 
the different from the ones in the MSSM. These processes may be used as the probe of 
the model in the future experiments. 

Branching ratio of the lepton flavor violating process £ a — > £p r ) is given as 



i=l a=± 



M}\ a _ / M„ 



2G F m 2 va \m 2 va J 2G F fi 2 v \ ^ 
xBr(£-^£-^u a ), (24) 



where m^ and M ia are the mass eigenvalues defined in eq. (IT4"|) and the statements below 
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Fig. 3 The left frame shows the contours for Br(/j, —> &y) (red solid lines) and Br(r — > ptrj) (blue 
dashed lines) under the conditions (f22|) imposed by the neutrino oscillation data. A thin black dashed 
line represents a line for m^„ = Mi. The right frame shows the mo dependence of each branching ratio. 
Thin black dotted lines represent the present experimental bounds for fi — > e'y and r — > fi'y, respectively. 
We use a GeV unit for Mi, mo and fi v . 



it respectively. The function F 2 (x) is defined as 

1 — Qx + 3x 2 + 2x 3 — 6a; 2 In x 
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(25) 



If we use the assumed flavor structure for the neutrino Yukawa couplings (117)1 . we find 
that 
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(26) 



Now we examine the parameter regions consistent with both the neutrino oscillation 
data and the lepton flavor violating constraints. By using the above formulas for the lepton 
flavor violating processes and the neutrino oscillation condition ( )22|) . the branching ratio 
of fi — > e7 and r — > /ry predicted by the model can be plotted in the (M 1; plane if 
M 3 , m and A are settled. It is useful to note that these bounds become more severe by 



12 



making A smaller generally since a smaller A requires larger values for \hi\ as found from 
eqs. (TT51) and (IT5|) . We fix these values to 

M 3 = 7.8 TeV, m = 0.48 TeV, A = 1.24 x 1(T 9 (27) 

as a typical example. The last one corresponds to X u Xd = 10~ 7 ' 7 and tan (3 = 16, for 
example. For a while, we consider the case ip± = ifi2 only. 

In the left frame of Fig. 3, we show the contours of these branching ratios. Red 
solid lines represent the contours of Br(fi — > cy) x 10 11 = 1.2 and 0.6 downward, and 
blue dotted ones are contours of Br(r — > fi'j) x 10 s = 0.1 and 0.04 rightward. The 
former one is independent of Mi as is clear from the expression in eq. (T26]) . F 2 (M|/m^ a ) 
becomes smaller for a larger M3 although the larger M3 makes the neutrino Yukawa 
coupling /13 larger through eq. ( 122ft . Thus, this branching ratio can easily satisfy the 
present experimental bounds by making M 3 large enough. We note that this feature 
is intimately related to the flavor structure (ITT)) which induces the tri-bimaximal MNS 
matrix. Br(r — > ^7) is found to show a different behavior, which is mainly controlled by 
the masses and the couplings of iV 12 . This figure shows that the present experimental 
bounds for these processes [35J are satisfied at the wide range of the m v _ > Mi regions, 
where = Mi is plotted by a thin black dashed line. 

In the right frame of Fig. 3, as in the left frame each branching ratio is plotted 
by a red solid line and a blue dotted line as a function of m for M% = 2.5 TeV and 
fi v = 4.2 TeV. The present bounds for them are also plotted by the thin black solid 
and dashed lines, respectively. The figure suggests that the constraint from the \i — V e^y 
requires mo ^ 0.45 TeV, although the r — > ^7 gives no constraint on tuq. From these 
figures, we find that the model can be easily consistent with both the neutrino oscillation 
data and the lepton flavor violating constraints for the natural values of parameters as 
long as they are suitably fixed. 

As a related subject, it is useful to note that the similar diagram to Fig. 2 contributes 
to the electric dipole moment of an electron (EDME) and the muon g — 2. Even if the 
neutrino Yukawa couplings are complex, the EDME is not induced since there is no 
mixing among AVs. New contributions to the muon g — 2 due to the similar diagram to 
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Fig. 2 are summarized as 

+ ~ 2 ((M 2 + \h 2 ?)F 2 + \h^F 2 f^f] 11 . (28) 

We estimate it for the allowed parameter sets obtained in the above analysis. The results 
seem to be smaller by three order of magnitude in comparison with <5a M = (30.2 ± 8.7) x 
1CT 10 , which is a discrepancy shown by the SM prediction and the value derived by the 
experiment [36]. This suggests that another origin is required for the explanation of this 
muon g — 2 discrepancy. 

2.4 Two dark matter candidate 

The model has two types of the DM candidate in general. One of them is the lightest 
neutralino x whose stability is guaranteed by the i?-parity as in case of the MSSM. The 
other one is the lightest neutral field with the odd parity of the new Z 2 symmetry, which 
is the remnant symmetry of the anomalous U(l)x- It corresponds to the lightest neutral 
state composed of the components of the chiral supermultiplets iVj, or rf^ d and 0. In the 
following study, we assume that the singlet fermion ip^ 1 (the fermionic component of N%) 
is the lightest one among these candidates. Since this Z 2 is not an exact symmetry but 
is weakly broken as shown in eq. (E]) by the anomaly effect, the latter candidate ip^ is 
not stable but it could have a long lifetime comparable to the age of the universe. The 
condition for this possibility is discussed in the next section. If this is the case, the DM 
relic abundance suggested by the WMAP [lj should be satisfied by these two contributions 
such as 

tt x h 2 + ttip Ni h 2 = 0.11. (29) 

The annihilation of two ipw^s is induced through the t- and w-channel r\ u exchange. 
Final states of such processes are composed of a pair of lepton and antilepton or a pair of 
slepton and antislepton. The latter final state gives new contributions compared with the 
non-supersymmetric long as Mi > m^(= mo) is satisfied for the slepton mass m^. 

If ipN 2 has a mass almost degenerate with the one of ipN 1 and this is the case here, we need 
to consider the coannihilation effect [37]. Final states of this annihilation are controlled 
by the flavor structure of the neutrino Yukawa couplings which is shown in eq. (flTj) . On 
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Fig. 4 Diagrams contributing to the i-channel annihilation of the two DM candidates. There are also 
u-channel diagrams corresponding to these. Although the hnal state // of the x annihilation is composed 
of both qq and i£, the annihilation contains £1 only. 

the other hand, the annihilation of two x' s occurs through various processes depending 
on its composition, which is determined by both the supersymmetry breaking scenario 
and the radiative symmetry breaking conditions (see Fig. 4). It has been studied in detail 
in the MSSM context [3]. Here we do not fix the scenario but just confine our study to 
the pure bino case by assuming that fi and the masses of Higgs doublet and gauginos 
are tuned to satisfy the required condition, for simplicity. In this case, the annihilation is 
expected to occur dominantly through t- and w-channel sfermion exchange] 13 ! We estimate 
the annihilation cross section av for these DM candidates by expanding it as av = a + bv 2 
in powers of their relative velocity v |38j . 

For the singlet fermion ip^ , we need to take account of the coannihilation effect with 
ipN 2 because of the assumption M\ ~ M^. In order to estimate the freeze-out temperature 
Tf of ^jvi including the coannihilation case, we follow the procedure given in [37J. We 
define a e s and g e ff as 



0"eff 



2 

ihh , r , 9N i gN 2 

i/eff y e ff 

g Nl +g N2 (l + 5) 3 / 2 e- x6 , 



5f' 2 e~ x6 + 



9jh 

9 a ^N 2 ll>N 2 

g c s 



1 + 8fe 



-2xS 



(30) 



where internal degrees of freedom of A7j are described by g^. and 5 = (M 2 — Mi) fM\. If 
we define a e s and 6 e ff by a c sv = a c s + b c sv 2 , the thermally averaged cross section can be 
written as (a e fiv) = a e ff + 6b cfi /x where x = Mi/T. Since 5 1 is supposed here, the 
second and third terms can bring the important contribution. Using these formulas, the 



13 If the lightest neutralino x i s heavier than Higgs scalars, we need to take account of the t-channel 
Higgsino exchange process. However, it is expected to be subdominant as long as the Higgsino is heavier 
than the sfermions, which is assumed throughout this analysis. 
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effective annihilation cross section is given by 

Y 4 M 2 

1* V M ii M t + m la m lb) 2 
96 7 r a ^ ± (M 1 2 + m2 a )2 (M 2 +m 2 b)2 ^ 

327r(/x 2 + M 2 / 3 2 ) 2 V 3 3(/x 2 + M 2 /? 2 ) V ' 1 J 

where (3 = y/l - (mo/M^ 2 . The mass eigenvalue m^ a is given in eq. ( fT4l) . K s and 1^, are 
defined by 

Y 4 = 2\h 1 \ 2 \h 2 \ 2 sin 2 ^! - <p 2 ), Y 4 = \h\ 4 + 2\h 1 \ 2 \h 2 \ 2 cos2(<^ - <p 2 ) + \h 2 \ 4 . (32) 

The first and second lines in eq. (l3Tj) represent the contributions with the lepton-antilepton 
final state. The third and forth lines come from the slepton-antislepton final state. We 
find that there can be s-wave contributions in the first and third lines if the neutrino 
Yukawa couplings have the phases such as (pi ^ ip 2 + ri7ro This happens since the 
coannihilation cross section <j$ N ^ N allows the s-wave contribution [7]. 

On the other hand, the s-wave contribution to the annihilation cross section of the 
bino-like % to // is expressed as [38] 



where z = mf/m x , (3' = (1 — z 2 ) 1 / 2 and Qj is the electric charge of /. In the sum- 
mation in eq. (133]) . / runs among quarks (cf = 3) and leptons (cf = 1) which satisfy 
m/ < m x . Although XX can annihilate to a pair of Higgs scalars through the Higgsino 
exchange, we neglect it as a subdominant process. If the bino mass M% is almost degen- 
erate with Higgsino mass fin, the coannihilation between x and a Higgsino-like neutralino 
x' brings important effects on the relic abundance of x- The relevant cross sections to 



14 The s-wave contributions are dominant in such cases. The expressions for pf,s f° r the corresponding 
p-wave contributions are shown in Appendix C for the completeness. 
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Fig. 5 The left frame shows the contours of 0^ Ni /i 2 + Q x h 2 = 0.11 for m = 0.46 TeV (a green dotted 
line), 0.48 TeV (a blue dotted line), and 0.50 TeV (a red solid line). Each figure in the right frames shows 
the ratio of each DM component to the total relic abundance h 2 /0.11 and fl x h 2 / 0.11, respectively. 
Each line corresponds to the same value of mo used in the left frame. In these figures ip\ — ip% = is 
assumed. 



this coannihilation are given by 



a xx ,v ~ 



°x'x' v - 



(g 1 cos2/3/<7sin j3)m w rrif^f (T 3L tan 8 ] w )m w mf(f 



[K r 4m^ + m ^ryv2 

(£f cot(3/2)(m x > + fi H )m f 
[{m 2 



m o + ft' 2 ™ 2 



x 

-1 2 



4m2)2 + m^r^] 1 /2 4(m 2 + p' 2 m 2 x ) 



(34) 



where £/ = cot (3 and Q = for / in the up sector and £/ = tan /3 for and also (f = 
for / in the down sector. The effective cross section including the coannihilation can be 
determined by the similar formulas to eq. fl30|) . The dominant contributions for it are 
expected to come from a channel with the tt final state if m x > m t is satisfied. If we focus 
our numerical study to the case with m x ~ m' x and fix both tan and m 2 H0 to suitable 
value^fl, these cross sections are determined by two free parameters m and m x . They 
should satisfy some required conditions. Since % should be lighter than the left-handed 
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In this study we assume a sufficiently large value for wv^o such, as ui^o 4?ri^., for simplicity. 
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sneutrinos, m x < m has to be satisfied. The new contribution to the \i — > e'y imposes 
the lower bound on m as found from Fig. 3. In relation to this constraint it is useful to 
remind that a larger M3 allows smaller values for mo- 

Now we examine the possibility to realize the required relic abundance by these two 
DM. If we follow the ordinary method given in [HTJEB], we can estimate the relic abundance 
by using the results for the effective annihilation cross section given above. Each relic 
abundance and Q x is given by the formulas 

n , 2 1.07 x 10 9 Xf _ 0.038g eS m p im DM (a C E + 6b cS /x f ) 

gj m p i(GeV)(a e ff + 3b e $/Xf) g* x/ 

where m v \ = 1.22 x 10 19 GeV and tudm is the mass of DM. g e n stands for the internal 
degrees of freedom of DM or the effective degrees of freedom in the coannihilation case. We 
can use ~ 100 as the relativistic degrees of freedom at the DM freeze-out temperature 
Tf(= todm/x/) for both ip^ and x- 

The results are shown in Fig. 5 for the parameter set used in the previous section. In 
the left frame, we plot the contours for the relic abundance (1291) in the (Mi,m x ) plane 
for niQ = 0.46 TeV (a green dotted line), 0.48 TeV (a blue dotted line), and 0.50 TeV (a 
red solid line). Since we are considering that \ is the DM lighter than Mi, m > m x 
should be satisfied. The allowed regions in the (Mi,m x ) plane are represented by the 
points on each contour which satisfy this condition. They are found to have the almost 
fixed values of m x for the larger values of M\. The reason can be found in the right 
frame of Fig. 5, where the relic abundance of each DM component is plotted for the same 
parameters as the left frame. The same lines are used as the ones for the corresponding 
contours in the left frame. Since Q^ N is almost constant at large Mi regions, the condition 
fl29|) can be satisfied only for restricted values of Q x . Since both DM components have 
the same order abundance, we can expect rather different DM phenomenology from the 
one component DM models. 

We need to note that the above solutions are obtained for rather large neutrino Yukawa 
couplings such as a/I^-iI 2 + I ^2 1 2 — 3.56 in case of m = 0.48 TeV, for example. Although 
these results on the neutrino masses and the DM are interesting enough for the model, 
such large Yukawa couplings are dangerous for the perturbatibity and the stability of the 
model [T7]. However, this point can be improved by considering the case (pi — (p% 7^ nir, 
which makes the s-wave contributions to the ^1 annihilation cross section possible. Since 
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Fig. 6 The left frame is the same one as the left frame in Fig. 5. However, ipi — ip2 — 7r/4 is assumed 
here. The right frame shows the phase dependence of Q^h 2 for each value of mo in the left frame. 

the s-wave contributions give much larger effects than the p-wave ones for the thermally 
averaged annihilation cross section (av) , the neutrino Yukawa couplings required to reduce 
the relic abundance of ^jvi can be smaller. In the left frame of Fig. 6, the contours 
corresponding to the ones in the left frame in Fig. 5 are plotted by assuming \h\\ = \hz\ 
and ipi — if2 = 7r/4. The parameters are fixed to the same values as in Fig. 5 except for 
A, which is taken here as A = 8.4 x 10~ 9 . Since the neutrino Yukawa couplings have much 
smaller values like a/I^-iI 2 + | 1 2 — 1.45 in case of mo = 0.48 TeV, the above mentioned 
tension is relaxed. In the right frame of Fig. 6, the phase dependence of Vl^^h 2 is shown 
for the same parameter setting. 

3 Probing two dark matter 

The present model has two DM components, that is, the meta stable lightest singlet 
fermion ip^ and the lightest neutralino \- They are expected to be observed by sev- 
eral kinds of experiments. The former one may be studied indirectly through the decay 
products such as charged particles and gamma rays. The anomaly reported in the cosmic 
rays by PAMELA and Fermi-LAT may be relevant to this decay. The latter one may 
be observed directly through the elastic scattering with nuclei in the same way as the 
ordinary lightest superparticle. However, the situation can be rather different from the 
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Fig. 7 Decay processes of V'jVi to x- A bulb represents the anomaly induced interaction CjBM p \e bj Ljfj u . 

MSSM since the DM relic abundance is composed of two components. In this section 
these subjects are briefly discussed and the detailed study will be presented elsewhere. 

3.1 Decay of the right-handed neutrino dark matter 

The relic DM in this model can be composed of the two components ip^i and x as shown 
in the previous part. However, one of the DM candidates ip^ is not stable since the Z 2 
symmetry which guarantees its stability is not exact. It should be reminded that this 
symmetry is the remnant symmetry left after the spontaneous breaking of the anomalous 
U(l)x- Green-Schwarz anomaly cancellation mechanism induces the Z 2 violating interac- 
tion and also the corresponding soft supersymmetry breaking term, which are shown in 
eqs. © and f flUj) . If ip Nl is heavier than x, this interaction brings the decay of ip^ to x 
through the diagrams shown in Fig. 7. 

We can estimate the lifetime of t/Vi due to the decay derived by this interaction. It 
can be expressed as 

/lTeV\ 5 / fi„ \4/ m \4 (\ TeV\ 2 ( e 2b > \ 2fi . , 

where we use hf x = 0(1) and Cj = 0(1). If Mi ^> m is satisfied, m should be replaced 
by Mi in eq. (|36l) . From this formula, we find that can have a sufficiently long lifetime 
compared with the age of universe as long as bj > 82 is satisfied. Thus, although the true 
stable DM is the lightest neutralino Xi we need to take account of the contribution of 
to the relic DM abundance in the universe as discussed in the previous part. 

On the other hand, depending on the scale of the Z 2 breaking CjM v \Be~ hi in eq. fflOl) . 
particles yielded in the decay of ip^ could bring additional contributions to the cosmic 
rays in the present universe. They could be detected as the anomaly in the expected 
flux of the cosmic rays through the various observation. In fact, if bj ~ 92 is satisfied, 
this anomaly induced superweak interaction causes a large enhancement factor of O(10 80 ) 
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Fig. 8 Elastic scattering diagrams of the bino-like neutralino \ with the quark which are relevant to the 
direct search of %• 



in eq. ( 13 6 p to realize a long lifetime of O(10 26 ) sec for ib^ J 16 l This lifetime is known to 
be suitable to explain the charged cosmic ray anomaly reported by PAMELA [12] and 
Fermi-LAT [13]. Since ip^ couples only with leptons and sleptons because of the U(l)x 
symmetry, this decay can yield a pair of a lepton and an antilepton, or a photon in 
addition to the lightest neutralino x- This feature is favored by the lepto-philic nature 
of the PAMELA observations. Moreover, the flavor structure of the neutrino Yukawa 
couplings (I17p can restrict the flavor of the final charged leptons to \x and tO Since this 
makes positrons and electrons produced by their decay soft, the model becomes favorable 
for the explanation of a plateau at high energy regions of e + + e~ spectrum found in the 
Fermi-LAT observations [39] . 

The heavier DM component has also a radiative decay mode to the lightest neu- 
tralino x- Its one-loop diagram is shown in Fig. 7. This decay associates a characteristic 
gamma which can be detected through the observation of the cosmic gamma rays. It is 
expected to appear as a line shape spectrum at the energy = (Mf — m'^)/2M 1 , which 
corresponds to the endpoint of the gamma-ray spectrum caused by the bremsstrahlung 
and the inverse Compton scatterings associated to the charged decay products of ip^- 
This could be a clear evidence of the model. 



3.2 Direct detection of the neutralino dark matter 

Direct detection of the DM is expected to clarify the nature of DM [3]. Several experiments 

such as CDMSII, XENON100 and XMASS to observe its elastic scattering with nuclei are 

16 It is interesting to note that this value of bj can be consistent with the value of TrX for gx = 0(1), 
which is required to realize the vacuum with the desired values of e± as discussed below eq. ([5]) and in 
Appendix A. 

17 In that case we need to impose c e = additionally in the anomaly induced interaction CjM p xe * Ljr/ U . 
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now under going or will start in the near future. In the study of DM models, it is crucial 
to address the discriminative features of the model, which are expected to be shown in 
these experiments. 

Since one DM component interacts with the leptons only and can not have in- 
teractions with nuclei at tree level, the scattering cross section with nuclei is heavily 
suppressed by the loop factor. Thus, it is difficult to detect it in these experiments. On 
the other hand, the neutralino DM \ can be scattered with nuclei at tree level since it has 
the same nature as an ordinary neutralino in the MSSM. One may consider that there is 
no distinction with the MSSM case. However, it should be noted that the detection rate 
R in this model can be different from the direct detection rate -Rmssm in the MSSM even 
if the DM x has the same mass and the same scattering cross section with a nucleon in 
two models. The detection rate is approximately defined by 



where rij is the number of % nuclei species in the detector and (<Ji X ) is the scattering cross 
section of the x an d % nuclei species averaged over the relative velocity between the x an d 
the detector. Since the present model has two DM components, the detection rate in this 
model is related to the MSSM one as 



for the fixed m x and m,Q. Therefore, in the present model the parameter regions to realize 
the same DM detection rate can be changed from the ones in the MSSM, although the 
interactions of the x with quarks are same as the MSSM neutralino. This might open a 
new possibility for the supersymmetry breaking parameters, which is considered not to 
be allowed in the MSSM case. 

In order to see this aspect we consider the bino-like x-, as an example. The dominant 
contributions to the elastic scattering of the bino-like % with nuclei come from the squark 
exchange shown in Fig. 8. In the MSSM case, the allowed parameter regions relevant 
to the detection rate have been extensively studied. We are interested here in possible 
changes of the detection rate from the MSSM and also the change of the allowed region 
of the supersymmetry breaking parameters. 

We suppose the values used in Fig. 5 for the parameters relevant to the neutrino mass 
generation. Eq. (138)) shows that the figure of R/Rmssm a s a function of m x is the same 




(37) 



R 




-Rmssm, 



(38) 
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Fig. 9 The Contours of R/Rmssm in the (m x ,mo) plane. Each line corresponds to R/Rmssm = 1 (a 
red solid line), 0.6 (a green dotted line), 0.5 (a blue dotted line) and 0.4 (a red dashed line). A thin black 
dashed line represents a line for tuq = m x . 

as the one in the right frames of Fig. 5 for the case <pi — <p2 = 0. If we considered the 
case with m = 0.48 TeV and Mi = 2.5 TeV as an example, Fig. 5 shows that m x should 
be 0.2 TeV or 0.45 TeV to satisfy the WMAP constraint and then R/Rmssm ^ 0.5 or 
0.47 for each m x value, respectively. These suppressed detection rates compared with 
the MSSM are realized for much smaller mo values than the ones in the MSSM. The 
values of (m x , m ) found through the direct search of the lightest neutralino are placed 
on the contours of Q x h 2 , which exist in the region where Q x h 2 < 0.11, mo > m x and 
tuq > 0.45 TeV. The last condition comes from the \i — Y ej constraint in Fig. 3. We note 
that the model could bring new possibilities for the detection rate of the neutralino DM 
and values of (m x , mo) relevant to the supersymmetry breaking even if we impose the DM 
relic abundance constraint (see Fig. 9). 

4 Summary 

We have studied a supersymmetric model with an anomalous U(l) gauge symmetry from 
a view point of the neutrino masses and the DM. The model considered in this paper 
may be recognized as a supersymmetric extension of the radiative seesaw model for the 
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neutrino masses proposed by Ma. The Froggatt-Nielsen mechanism based on the sponta- 
neously broken anomalous U(l) symmetry generates the hierarchical structure of Yukawa 
couplings of the quarks and the charged leptons. Moreover, it can also explains the 
hierarchical couplings and mass scales required for the radiative generation of the neu- 
trino masses. If we assume a flavor structure for the neutrino Yukawa couplings, the 
tri-bimaximal mixing is automatically induced. 

The model has two DM components. One is stable and the other is the decaying 
DM. Their stability is guaranteed by two Z 2 symmetries, one of which is the ordinary 
R parity. Since one of these discrete symmetries is assumed to be anomalous and then 
the Z 2 violating interaction is generated nonperturbatively, the instability of one DM 
component is caused. This phenomenon happens since this Z 2 symmetry is embedded into 
the anomalous U(l) gauge symmetry. Since the Z 2 violating nonpertubative interaction 
is extremely weak, the huge suppression factor for the decay width of this DM can be 
derived. As a result, its lifetime becomes longer than the age of the universe. These DM 
particles can be detected through both direct and indirect searches. Thus, the model may 
be checked through these experiments. In particular, the recent and future data coming 
from the cosmic ray observations might clarify the relation between this model and the 
physics beyond the SM. 
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Appendix A 



In this appendix we give two examples of the U(l)x charge assignment which induce the 
favorable effective parameters for the generation of the mass eigenvalues and mixings of 
the quarks and the leptons including the neutrinos. If we define e± as e± = (E±)/M p i, 
eq. (J7J) gives the expression for each parameter in the superpotential W such as 



h U _ V n Qi +n Uj+ n Hu 

'Hi - y%j L - ■ 



\i H = X H M pi e_ 



(39) 



for the ones belonging to the MSSM. Other ones are also expressed as 



\ ,,.V £ n L i + n N j +n r , u +l 

n Vu+ n V d + 1 



K = vS e 



n Vu+ n <t,+ n H d + 1 n Vd +n <t> +n H u + 1 

= Vrtu , Ad = y Vd e_ 



n N i + n N j + 1 



Mij = y NlNj M pl e + 



^ = y^M pl e 2 y +1 . (40) 



In these formulas, it may be natural to suppose that the original coupling constants in 



the nonrenormalizable interactions, that is, y¥ Xh and so on, have values of 0(1). 

We assume that XL and S + obtain the VEVs defined by £_ ~ 1CT 1 and e + ~ 
respectively. The U(l)x charge is assigned to each field as follows, 



10" 



n Q = (6, 5, 3), n v = (6, 5, 3), n D (4, 3, 3), 
example (i) n L = (9,9,9), ^ = (1,-1,-3), n Hu = n Hd = -6, 

n N = (-5, -5, -5), n Vu = n Vd = -5, = 7, n + = 3. 

In this case, it is easily found that the mass matrices for the up- and down-type quarks 
and the charged leptons take the following form: 



/ 



Mi 



u 



e 5 e 



3 \ 



e 5 si e 2 



\e 3 _ el 1 / 
/ e 4 \ 



el 

V o o i ) 



(H° u ), 



M, 



D 



el el 



e. \ 
1 

1 J 



(H° d ), 



(41) 



where these Mf are defined as ipuMfipL- In the charged lepton mass matrix, the off- 
diagonal couplings are supposed to satisfy y? = (for i ^ j). From these mass matrices 
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we obtain both the mass eigenvalues and the CKM matrix in the quark sector as 
m u : m c : m t — e Q _ : si : 1, rrid : m s : mj = st_ : e 2 _ : 1, 
K s ~£-, K6~ei, V cb ~5 2 . (42) 

The charged lepton mass eigenvalues satisfy 

m e : m M : m T = e A _ : e 2 _ : 1. (43) 

These can give qualitatively good results as long as takes a value of Cabibbo mixing 
angle 0.22. The effective neutrino Yukawa couplings have no hierarchical structure = 
yfj = 0(1) as supposed in the text. Since the effective \xh term is generated as M v \e A +) it 
can take an appropriate value for the electroweak symmetry breaking. 
Other effective parameters are estimated as 

A u ~A d = 0(e+), M t = 0{y NiNt M pl e 3 + ), //„ = 0{y v M pl e 3 + ), ^ = 0(M pl ei 5 ), (44) 

where we suppose that the off-diagonal couplings y at 4 at,- are zero. These parameters are 
intimately related to the neutrino mass generation in the present model. The values of 
these parameters used in the text can be realized if yjVjJVj and y v have suppressed values 
of 0(1CT 2 ). 

In the above example, we implicitly assume that each parameter is determined by 
either e_ or e + only but is not determined by both of them. In the next example, we 
consider that some of them are determined by both e_ and e+. We assume that £_ and 
S + obtain the VEVs defined by e_ ~ 1CT 1 and e + ~ 4 x 10~ 4 again. The U(l)x charge 
is assigned to each field as follows, 

n Q = (6, 5, 3), n v = (3, 2, 0), n D = (5, 4, 4), 
example (ii) n L = (7, 7, 7), n E = (4, 2, 0), n Hu = -3, n Hd = -7, (45) 
n N = (-6, -6, -6), = -2, n m = -9, = 6, n+ = 3. 

Using these charges, we can easily find what kind of factor determined by E- and e + 
should appear as the lowest order one for each term in W. The mass matrices of quarks 
and charged leptons and also the neutrino Yukawa couplings hfj show the same features 
as the ones in the previous example. 

The remaining effective parameters are estimated as 

Hh = (M pl eXe 2 _ ) , X u = O (e+e„ ) , X d = O (e 2 + e- ) , 

Mi — 0(M pl £^£_), N = 0{M vX e%e 2 _), ^ = 0{M pl e 13 ). (46) 
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These also result in the favorable values for the parameters relevant to neutrino mass 
generation and also the electroweak symmetry breaking. However, it is useful to note a 
following point. Since X u ^> X d is satisfied in this case, the formula for the dominant 
contribution to the neutrino masses can be changed. In fact, the term with X\ in eq. (14T)) 
could cause the similar order contribution as eq. (|T5|) or larger contribution even if there 
is a suppression due to the large fi^. 

Appendix B 

In this Appendix we give the formulas for the contributions to the neutrino masses due 
to the diagrams (b) and (c) in Fig. 1. These diagrams have the scalar component and 
the fermion component of the chiral superfield <fi as an internal line, respectively. The 
expressions for their contributions can be derived by taking account of eq. (I14p . 
The neutrino mass matrix induced by the diagram (b) can be expressed as 

i a a \ ST^ h ai hpiMi(A + a^) 2 a 

[M v u -2^2^ (4^2 

a=± i=l V ' 

X 2 v 2 ( 

i^I(m v+ ,m v+ , Mi, m^a) + I(m v _,m v _, Mi, m^ a ) - 2I(m v+ , m v _, Mi, m^J 

X 2 v 2 ( \ 
_l — |_^j( m??+)mj3+) Mj,m^ a ) + I{m ri -,m ri -,M i ,Tn^ a ) + 2I(m rj+ ,m ri _,M i ,m <j}a )J 

. X u X d v u v d a ( \ 
H t [I \m v+ ,m ri+ , Mi, m^) - I(m v -, m v ^, M h m^ a ) J 

where the function / is defined as 
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On the other hand, the neutrino mass matrix induced by the diagram (c) is estimated 
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as 



(M v )a/3 = n( u . ? 9 h ai hpii(^ hjj, M ia , fif) 

^ W a=± i=l 



\2 



(4tt) 

v ' o=± i=l 



y y > shgjhpjJtfi,,, fif), (49) 



where J is defined as 



J(m a ,m b ,m c ) = -■ — ^ — ^ + — ^ a b eZ- lnm^ 

(m^ — m A a )\m A c — m A a ) 2{m^ — m z a ) z [m z c — m a ) 
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^ / o Own <)\ o 1 



2(m^ — m\){m 2 a — m\) 2 2{m 2 — m 2 )(ml — m 2 ) 2 

2 



2(m2 - m 2)( m 2 - m 2)' ( 50 ) 
As long as /x^ 3> fi v , Mi, m , A and A u ~ are satisfied, these contributions are found 
to be subdominant in comparison with the one induced by the diagram (a) except for the 
ones proportional to f£ in eq. ( 14"T|) . These cause non- negligible contributions to eq. (jl5|) . 
which are estimated as 



Eh a ihpiMi 
X 2 v 2 



4 



I(m v +,m v+ ,Mi) + I(m v -,m v -,Mi) — 2I(m v+ , m v -,Mi) 



+ ^^(l(m v+ ,m v+ , MA + I(m v _,m v _, Mi) + 2I(m v+ ,m v _, M^ 
where the function I is defined as 

~. , m 2 him 2 . ml lnm? m 2 \nm 2 

I(m a ,m b ,m c ) = — ° a + — \ b ^ + 



(51) 



[m 2 - m 2 a ) (m 2 c - m 2 a ) (m 2 - m 2 ) (m 2 a - m 2 ) (m 2 - m 2 c )(m 2 a - m 2 ) ' 

m 2 — m 2 + m 2 \n(m 2 /m 2 ) , . 

I(m a ,m a ,m c ) = -2 c c \ c/ a> . (52) 

\m A a — m^y 

This could be the same order contribution to the neutrino mass matrix as the one given in 
eq.( TT5]) . However, it is useful to note that these contributions do not change the structure 
of the MNS matrix but they only change the values of Aj somewhat. 

Appendix C 

The p-wave contributions in case of (pi ^ ipi + nix are given by 

1 M 2 (3jg + 5mgJ M 2 (3M 2 + 5m^) Mf 

PF ~ 3 " 6(M? + m 2 va ) 2 6(M 2 + m 2 b ) 2 + 3(M 2 + m 2 a )(M 2 + m\ h ) ' 

1 < 5 o2 , M^ 2 (^ + 4MOT 

2 + 12^ + 3(/i2 + M 1 2 /32)2 ■ ^ 
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